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ABSTRACT 

Astrometric monitoring of the Sirius binary system over the past century has yielded several predictions for 
an unseen third system component, the most recent one suggesting a <50Mjup object in a ~6. 3-year orbit 
around Sirius A. Here we present two epochs of high-contrast imaging observations performed with Subaru 
IRCS and AO 188 in the 4.05 /im narrow-band Br ol filter. These data surpass previous observations by an order 
of magnitude in detectable companion mass, allowing us to probe the relevant separation range down to the 
planetary mass regime (6-12 Mjup at \" , 2-4Mjup at 2" , and 1.6Mjup beyond We complement these 
data with one epoch of M-band observations from MMT/AO Clio, which reach comparable performance. No 
dataset reveals any companion candidates above the 5 a level, allowing us to refute the existence of Sirius C 
as suggested by the previous astrometric analysis. Furthermore, our Br a photometry of Sirius B confirms the 
lack of an infrared excess beyond the white dwarf's blackbody spectrum. 

Subject headings: planetary systems — techniques: high angular resolution — stars: individual (Sirius A) — 
stars: individual (Sirius B) — white dwarfs 



1. INTRODUCTION 

One and a half centu ries ago, Sir ius was found to host a 
faint binary companion (lBondlll862l) . The orbital motion of 
this pair has been monitored ever since, leading to a num- 
ber of publications that claimed to find periodic perturbations 
indicative of the presence of an unseen third system compo- 
nent, Sirius C. The most recent analysis predicted a substellar 
companio n in a ^6. 3 -ye ar circumstellar orbit around Sirius A 
(feenest & Duven3ll995 . and references therein). While the 
amplitude of the purported astrometric signal, 56 mas, would 
suggest a companion mass of 72 Mjup, these authors imposed 
an upper limit of < 50 Mjup on the basis of system stability 
considerations. Since no measure of confidence is given, we 
assume a conservative lower limit of half the measured am- 
plitude, 28 mas, on the basis of their plotted periodograms, 
resulting in a minimal mass of 36 Mjup. 

This would place Sirius C in the so-called brown dwarf 
desert, the range of orbital parameter space around stars 
in which bro wn dwarfs are found to be s carce (e.g. 
iMarcy & B utler '20001: iGrether & Lineweav er 2006). Numeri- 
cal simulations of the formation and evolution of brown dwarf 
companions reproduce this scarcity, regardless of the forma- 
tion process assumed (e.g. as part of the star formation pro- 
cess (Bate & Bonnell 2005), planet formation by core ac- 
cretion (Mordasini et al. 2009) or gravitational fragmentation 
iStamatellos & Whitworth 2009)), thus these rare objects im- 
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pose important constraints on theory. 

However, Sirius' extreme brightness - both apparent and 
absolute - had long thwarted attempts to verify these claims 
through direct imaging. Kuchner & Brown (2000) estab- 
lished first constraints on substellar companions at separa- 
tions of V!5-3^^ with space -based coronography at 1 . 02 /im 
on HST NICMOS, whereas iBonnet-Bidaud & PantinI (120081) 
used ground-based observations assisted by adaptive optics 
on ESO ADONIS to achieve similar constraints in the range 
of 3" -10". These limits left most of the parameter space 
in which Sirius C was expected unexplored. In this work, 
we present new results from 4.05 /im observations on Sub- 
aru IRCS as well as from M-band observations on MMT/AO 
Clio, which improve the companion mass constraints by an 
order of magnitude and extend the coverage down to an inner 
working angle of O'H. With a baseline of 4.3 years among the 
observations, the chance that a companion be missed in all 
datasets due to adverse geometric conditions is slim. 

Sirius A is an AlV-type star at a distance of 2.64 pc, a mass 
of 2.02 Mq and an age of 225-250 Myr, whereas Sirius B is a 
white dwarf of 0.98 Mq with a cooling age of 124 ± 10 Myr 
orbiting Sirius A with a 50-year period (e.g. Lie bert et al.l 
l2005h . The combination of youth and extreme proximity (cf. 
a medi an target distance of 22 pc in the Gemini Deep planet 
survey. lEafreniere et al.ll2007bl ) allow us to explore unusually 
small orbital radii, with planetary-mass detection limits down 
to separations of 2.5 AU in projection. 

2. OBSERVATIONS 

The Sir ius system was observe d with the Subaru IRCS in- 
strument (iKobayashi et al.l [20001) on Ja nuary 20, 2011, us - 
ing the AO 188 adaptive optics system (|Hayano et al.ll201Q|) . 
The Br a narrow-band filter at 4.05 /im was chosen, since this 
band has been shown both theoretically and observationally 
to provide the best combination of PSF quality and brightness 
contrast for planet detection in th e speckle-dominated regime 
(iJanson et al.ll2008L I2009L 120101) . The dataset consists of 56 
frames, each comprising 240 co-adds of 0.1 s exposures, for 
a total integration time of 22.4 min. Although IRCS does not 
offer Lyot coronography, a mirror plate with several circular 
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Fig. 1. — (a-c) Results of the January 2011 ADI observations of the Sirius system in the Bra band with Subaru IRCS. (a) Signal-to-noise map of the inner 
r^lO" calculated in concentric annuli after convolution with a 4px diameter circular aperture. The stretch is [— 5cr, 5cr]. The location of Sirius A is 

marked with a white plus sign. The short red lines indicate the position angles where residuals of the spider diffraction pattern slightly elevate the noise level. 
The yellow circle marks the area of increased noise generated by one of the unused mirror holes in the field of view, (b) Full field of view of the ADI intensity 
image. Sirius B is visible to the East of Sirius A. (c) Non-differential image of Sirius B. Unsharp masking on the scale of 15 pixels (3 X/D) has been applied to 
reduce background structures. The point source just outside Sirius B's first Airy ring is a ghost caused by a filter. The cutout shows Sirius B's unsaturated PSF 
core, (d-f) The corresponding images for the March 2011 Subaru IRCS data. 
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holes is available to reject light from the primary star. We used 
a hole with a projected radius of 0V3 to avoid excessive satu- 
ration; nevertheless, the shoulders of Sirius A's point-spread 
function (PSF) locally saturate out to a radius of ~0V6. The 
plate scale was 20 mas per pixel. 

Follow-up observations were taken with Subaru IRCS on 
March 12, 2011. Using the same observing strategy, 83 
frames were taken for 33.2 min of total integration time. 

Furthermore, we make use of M-band imaging data 
taken on D ecember 5, 2006, with the Clio 3-5 fim imager 
(iFreed et al. 2004) in conjunction with the adaptive secondary 
mirror on the 6.5 m MMTO telescope OBrusa et al.ll20 04). The 
star was nodded 5'! 5 along the long axis of the detector after 
five images were taken. Each of the 322 images consists of 50 
co-added exposures of 209.1 ms length, for a total co-added 
duration of 56.1 min. To avoid variations in the pattern of il- 
lumination on the Clio detector, the instrument is fixed in ori- 
entation with respect to the telescope, resulting in total field 
rotation of 24.6° for our Sirius data. Conditions were photo- 
metric, and the seeing 0V5-0V7 (as seen with the video rate 
optical acquisition camera) at an air mass of 1.50-1.64. The 
plate scale was 49 mas per pixel. 

All observations were taken in pupil- stabilized image ori- 
entati on to enable the an gular differential imaging technique 
(ADI. lMarois et al.ll2006 ). We employed the LOCI algorithm 
(Local ly Optimized Combination of Images. ILafreniere et al.l 
12007 a l) in order to search for faint point sources in Sirius 
A's speckle halo. This form of ADI is the most power- 
ful high-contrast imaging method currently available, as ev- 
iden ced by recent direct detections of substellar companions 
(e.g. iThalmann et a l. 2009; Marois et al. 2010) and e ven cir - 
cumstellar disks (Ffhalmann et al..2010: .Buenzli et ani2010l) . 
We employed the LOCI parameters descri bed as optimal for 
the test dataset in 'Lafrenier e et al.l (|2007al) . with a frame se- 
lection criterion of 0.5FWHM to avoid self- subtraction of 
companion signals. Given the excellent resulting image qual- 
ity, we refrained from further parameter fine-tuning. In all 
datasets, the inner working angle of 0^7 is due to detector 
nonlinearity and saturation rather than insufficient field rota- 
tion. 

3. RESULTS 
3.1. Subaru IRCS Br a Data 

The final Br a images after ADI reduction are presented in 
Figures [T^-f . No obvious point- sources around Sirius A are 
visible. To confirm this numerically, we first convolve the 
image with a sampling aperture of 4 pixel diameter, and then 
calculate the signal-to-noise map (S/N) by dividing the pixel 
values in concentric annuli around the star by their standard 
deviation. As expected, we find no signal above the 5 a level, 
discounting the locally elevated noise in the shadow of the 
unused mirror holes. 

We do find a conspicuous signal in the immediate vicinity 
of Sirius B in the January 2011 data (Figure [TJ;). However, 
closer inspection of the signal's behavior in the time-series 
of images reveals it to be locked to the position angle of the 
pupil rather than the field, and identifies it as a ghost. No such 
signal is seen in the March 2011 data. Apart from this, the 
environs of Sirius B are background-limited outside the first 
Airy ring. 

We use the unsaturated PSF of Sirius B in the Br a images 
for photometric calibration. Although no prior Br a photom- 
etry of Sirius B exists, its overall brightness in the near- to 




Fig. 2.— Signal-to-noise map of the inner 4^^x4^^of the MMT/AO Clio 
M-band data. The stretch is [— 5 cr, 5 cr]. No significant signal is found. 



mid-infrared c an be predicted by a 25, 193 ± 37 K blackbody 
(iBarsto w et al." "2005; "Skemer & Closel 1201 Ih . Model spec- 
tra by iLejeune et al. (1997) for a range of log^ values up 
to 5 cm/s^ for Teff = 25,000 K indicate that while absorption 
lines are present in the mid-infrared, they exhibit rather small 
equivalent width. For the Br a absorption line we estimate an 
equivalent width of 0.7 nm, and hence conclude that it does 
not present a significant deviation from a blackbody. 

As an independent confirmation, we calibrate Sirius B's 
brightness using a brief observation of the A IV- type star 
HD 40138 as a reference. The result, 9.27 ± 0.17 mag, is con- 
sistent with our expected value of 9.17 ± 0.10 mag, which is 
based on the J-band photometry by Bonnet-Bidaud & Pantin| 
(I2008h and stellar spectra by CastelH & Kurucz (2003). 
This provides further evidence that Sirius B does not have 
an infrared excess fro m circumstellar dust as proposed by 
iBonnet-Bidaud & Pantin (2()08). 

3.2. MMT/AO Clio M-band data 

Individual images from Clio are beam- subtracted, and a 
sigma clipping routine is used to remove hot pixels deviat- 
ing by more than 3 a from a 5 x 5 pixel region centered on 
that pixel. The location of the central star is estimated by 
smoothing the images with a Gaussian kernel with FWHM 10 
pixels and estimating the location of the resultant peak. All 
the science images a re then passed i nto a LOCI data reduc- 
tion pipeline (La freniere et al.ll2007 al). The LOCI processing 
carries out reduction in 4 pixel (0^.'2) wide rings starting from 
an inner radius of 4 pixels (0'.'2). Each optimization section 
extends 20 pixels beyond the outer radius of the current ring, 
ensuring that the optimization section contains the equivalent 
area of a minimum of 150 PSF cores. After the LOCI back- 
ground subtraction is performed, the resultant images are ro- 
tated into the frame of the sky before being averaged together 
to give a final output image. 

Again, no significant signal above the 5 a threshold is 
found. The S/N map is shown in Figured 
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Fig. 3. — Limits on faint companions to Sirius A. (a) Effective 5cr de- 
tectable contrast to Sirius A in the January 2011 Bro; data, (solid red), the 
M-band data (dashed green), and the March 2011 Bro; data (dash-dot-dot- 
dotted magenta) as a function of separation, corrected for partial subtraction, 
(b) The companion masses corresponding to those contrasts, derived accord- 
ing to an updated version of the method of Janson et al. ( 2008) from the COND 
evolutionary model ( Allard et al. 2001; Baraffe et al. 2003). For comparison, 
the constraints established by Kuchner & Brown (2000, dash-dotted orange) 
and Bonnet-Bidaud & Pantin ( 2008, dotted blue) are overplotted. The gray 
shaded area indicates the parameter space of a <50Mjup companion in a 
6. 3 -year orbit (semi-major axis a = V/6). 



3.3. Contrast and detectable companion mass 

The 5 cr contrast curves established by the three datasets is 
presented in Figure [3] The partial self- subtraction of point 
sources during the LOCI reduction is estimated by inserting 
artificial 10 a point sources into the raw data across the en- 
tire usable range of separations from the star, and retrieving 
their surviving flux after subjecting them to the LOCI pipeline 
(radial filtering and optimized background subtraction). The 
plotted curves are corrected for this flux loss and therefore 
represent effective detectable contrast. The process is de- 
scribed in detail in lLafreniere et al.l (|2Q07a. Section 4.30: 

The COND -based evolutionary models (Allard et^ l2QQll: 
lBaraffeetal.ll200l are used to convert the M-band bright- 
ness contrast into detectable mass. For the B r a data, we 
adopt an updated version of t he method used in I Janson et aH 
(l2QQ8h . COND model spectra (lAflard et al.1 '2001') are used to 
re-calculate ph otometric prediction s from COND -based evolu- 
tionary models (iBaraffe et al. 2003) from L' band to the IRCS 
Bra band, for masses of 1-20 Mjup and the age of the Sirius 
system (250 Myr). The measured brightness contrast as func- 
tion of angular separation is then translated into detectable 
mass, using these values. 

The Bra data are sensitive to masses in the planetary 
regime down to the inner working angle of around Sirius 
A (6-10 Mjup at V\ 2-3 Mj^p at 2"), and reach 1.6 Mj^p in the 
thermal background beyond A" , which includes the immedi- 
ate surroundings of Sirius B down to 0'.'2. The M-band data 
achieve comparable performance in a more restricted field of 
view «3'.'5). 



4. DISCUSSION 

Our high sensitivity to substellar companions around 
Sirius A allows the first thorough test of the hypothesis of a 
<50 Mjup body, '' Sirius C", in a 6. 3 -year o rbit around Sirius 
A as proposed by lBenest & Duveni (Il995|) . As these authors 
point out, though, a single negative detection does not prove 
the non-existence of such an object. Even though the semi- 
major axis of such an orbit is a = 4.3 AU, corresponding 
to 1'.'6 at the distance of 2.64 pc, both projection effects 
and eccentricity can leave the companion at an apparent 
separation below the inner working angle of O'.'T at the time 
of observation. 

To quantify the likelihood of a false negative result, we sys- 
tematically calculate the projected orbital ellipses for 6. 3 -year 
orbits with eccentricities e = {0,0.04,.. .,0.96} and argu- 
ments of periastron uj = {0°, 10°, . . . , 350°}. Since we make 
use of radial detectable mass curves (Fig. [3]), we are insensi- 
tive to the orbit's overall position angle on sky; thus no vari- 
ation of the longitude of the ascending node ft is necessary. 
As for the inclination, we explore two scenarios: The copla- 
nar case, where Sirius C shares the orbital plane of the Sirius 
AB system with an inclination of 136.6°, and the free case, 
where all system orientations on the unit sphere are consid- 
ered equally likely, resulting in a statistical weight of sin i for 
each inclination i = {0°,1°,...,90°}. The mass of Sirius C 
is sampled as m = {0.5, 1, . . . , 20} Mjup. 

On each ellipse, 1000 possible companion positions at the 
epoch of our M-band observations are chosen, with uniform 
spacing in eccentric anomaly. We then calculate the corre- 
sponding location at the epoch of the other two datasets. The 
detectable mass curves are evaluated to determine whether or 
not the companion would have been detected at > 5 a at least 
once in that configuration. For each ellipse, the detection like- 
lihood is averaged over all companion positions, whereby a 
statistical weight is assigned to each position to account for 
the non-uniform evolution of eccentric anomaly with time. 
The weight is proportional to the time that the companion 
spends at each position over the course of an orbit, approx- 
imated by half the distance to the two adjacent positions in 
the sequence divided by the local orbital velocity. Finally, the 
detection likelihoods are averaged over all i and uj to yield the 
detection completeness for each combination of eccentricity 
e and companion mass m. The results of this analysis are 
plotted in Figure (H 

We find that coplanar companions down to 6 Mjup can be 
excluded at the 5 a level at 100% completeness. The com- 
pleteness remains above 50% down to 2-3.5 Mjup, depending 
on eccentricity. If no constraints are imposed on the inclina- 
tion, edge-on orbits emerge that can hide Sirius C behind Sir- 
ius As glare. As a result, the completeness values for masses 
above 12 Mjup drop to 97-99% for certain eccentricity ranges, 
and down to 90% for masses of 5-7 Mjup. 

The astrometric signal reported by iBenest & DuventI 
(|1995|) . 56 mas, implies a companion mass of 72Mjup. The 
authors furthermore impose an upper mass limit of <50 Mjup 
on the basis of system stability considerations. While no 
lower mass limit is given, we derive a conservative estimate of 
half that value, 36 Mjup, from their published periodograms. 
Given our much lower detection limits, our three combined 
epochs of high-contrast imaging can therefore decisively re- 
ject their Sirius C hypothesis. 

Although Benest & Duvent (|1995[) provide no error esti- 
mation, we consider a false alarm the most likely explana- 
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Fig. 4. — Probability of at least one 5 a detection of a companion with a 
6. 3 -year orbit around Sirius A in our three datasets, as a function of compan- 
ion mass m and eccentricity e. (a) The case of coplanar orbits with the Sirius 
AB system (inclination i = 136.6°). All companions above 6Mjup can be 
rejected regardless of orbital phase (100% completeness). The completeness 
remains above 50% down to 2-3.5 Mjup depending on eccentricity, (b) If the 
inclination is unconstrained (p{i) cx: sin i), edge-on orbits become possible, 
introducing a chance of a few percent that Sirius C might elude detection in 
all observations. Nevertheless, our three datase ts provide decisive evide nce 
against the existence of Sirius C as proposed bv lBenest & DuvenB (119951) . 
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Fig. 5. — Probability of at least one 5 a detection of a companion around 
Sirius A in our three datasets, as a function of companion mass m and semi- 
major axis a. A flat distribution in eccentricities is assumed, (a) The case of 
coplanar orbits with the Sirius AB system (inclination i = 136.6°). (b) The 
case of unconstrained inclination (p(i) oc sin i). The false negative probabil- 
ity^ large for the stable planet orbit regime at a < 17 (iHolman & WieserS 

tion for their results. Precision astrometry is known to suffer 
from a multitude of systematic errors (e.g. subtle changes in 
pixel scale and orientation, and differential atmospheric re- 
fraction which depends on airmass, parallactic angle, and am- 
bient conditions), and has led to a series of spurious detec- 
tions in the past (e.g. Pravdo & Shaklan (2009) vs. Bean et al. 
12010) and Lazorenko et al. (»2011i) : Ivan de KampI (119691) vs. 
[Gatewood & Eichhom (1973^). 

Leaving astrometric predictions aside, we can also explore 
the parameter space for other semi-major axes. Figure |5] 
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shows the completeness as a function of semi-major axis 
a = {0.25, 0.50, 10.0} AU and companion mass m, 
assuming a flat distribution in eccentricity. Although our 
data are sensitive to planets down to an inner working 
angle of ff!7 and down to 1.6Mjup at large separations, the 
completeness values drop quickly at shorter separations (e.g. 
50% for a lOMjup object at a = 2AU). This coincides with 
the domain of lon g -term stable planet orbits predicted by 
iHolman & Wiegerd (Il999h . with a critical semi-major axis 
ttc = 2.17 AU. Therefore, plenty of parameter space remains 
for unseen planets around Sirius A. The upcoming next 
generation o f high-contrast instruments, such as SPHERE 
(iBeuzit et"an[20 10). will offer smaller inner working angles 
and better contrast performance, and thus stand a good 
chance to detect such planets. In particular, due to its extreme 
proximity and brightness, Sirius is the third most promising 
target (after a Centauri A and B) for the direct detection 
of exoplanets in refl ected Hght with the SPHERE ZIMPOL 
imaging polarimeter (IThalmann et al.i r2008). 

One thing to keep in mind is the fact that Sirius B was orig- 
inally a ~5 Mq progenitor star that expanded into a super- 
giant ~125 Myr ago, with potentially dramatic consequences 
for the system architecture (Liebert et al. 2005). Accretion 
of ejected material from Sirius B may have caused planets 
around Sirius A to gain mass and heat, migrate, o r form 
in the first place as second-generation planets (e.g. iPeretsI 
"2010"). Since these processes would leave the planets hotter 
and brighter than their unperturbed 250 Myr-old counterparts, 
our detectable planet mass curves in Figure |3]are conservative 
for such objects. 

5. CONCLUSIONS 

We present three high-contrast imaging datasets of the Sir- 
ius system, all of which reach detection performances in the 
planetary regime (6-12 Mj^p at V\ 2-4Mjup at 2", 1.6Mjup 
beyond 4^0- This constitutes an improvement of an order of 
magnitude in detectable planet mass. Taken up to 4.3 years 
apart, the observations allow us to refute the existence of a 
substellar companion with a mass of <50Mjup in a 6.3-year 
orbit as predicte d from astrome try measurements of the Sir- 
ius AB system (iBenest & Duventlll995l) . For a companion 
mass above 12Mjup, the chances of a triple false negative at 
a 5 cr threshold are 0-4%, depending on eccentricity. For the 
special case of coplanar orbits, the probability is 0% down to 
6 Mjup). However, we note that our observations leave open 
the possibility for Jupiter- and Neptune-sized planets around 
Sirius A, especially at short angular separations. 

Furthermore, we confirm the absence dust around Sirius B 
by the lack of an infrared excess at 4.05 /im within our preci- 
sion of 0.17 mag. 

We thank David Lafreniere for generously having provided 
us with the source code for his LOCI algorithm. We are grate- 
ful for the privilege of observing from Mauna Kea, which 
holds great cultural significance for the Hawai'ian indigenous 
community. 

Facilities: Subaru (IRCS, A0188), MMT/AO (Clio). 
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